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Abstract 

In this talk I will introduces two spaces: the first space is the usual n-dimensional 
vector space with the unusual feature that n is non-integer, the second space is composed 
by the linear matrices acting on the previous space (physicists are particularly interested 
to study the limit n going to zero). These two spaces are not known to most of the 
mathematicians, but they are widely used by physicists. It is possible that, by extending 
the notion of space, they can become well defined mathematical objects. 

1 Introduction 

In last thirty years physicists have commonly used spaces having non-integer dimensions. 
These spaces do not have a clear mathematical meaning and it not evident in which sense 
the name space is not abusive. On the other end it would be very interesting if one could 
properly define these objects, that are commonly used by physicists. It is possible this task 
may be achieved by generalizing the notion of space. 

There are many examples of the previous construction. Here we will concentrate our 
attention on two different cases: 

• The usual space R n , when n is no more a positive integer, but it is a real number 

• The space of matrices n x n, for n real. In this case we are mainly interested to study 
in the limit n->0 |. 

The strategy followed by physicists in denning these strange objects consists in picking 
some properties of bona fine spaces for positive integer n and in making the analytic contin- 
uation to non integer n, without paying to much attention to the definition of the objects 
they are using. 

A clarification of the mathematical structures involved would be extremely interesting 
because these spaces are extremely useful and their introduction has a strong heuristic value: 
it allowed us to obtain many results that later on have been proved using conventional 
techniques. It is unlikely that this task can be done without generalizing the notion of space 
in an appropriate way. 
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2 The case of R n with non-integer n 

One of the physical motivations for introducing spaces with non integer dimension is the 
following. There are some physical systems whose properties cannot be computed directly 
in the case of the three dimensional space R 3 , however one can formally compute these 
properties in the space R n (that has dimension n) for n near to some specific value (e.g. 
n = 4-e Hi. 

The procedure used by physicists for introducing non-integer dimensional spaces is very 
simple. The properties we would like to compute can be expressed in terms of some integrals. 
In some cases the integrals can be computed for any dimension and their value can be 
extended to a meromorphic function of the dimensions [[|. 

Let us consider a very simple example: the volume of the <i-dimensional unitary sphere 
is given by 

vii) ~ nff^T) • « 

The previous formula for V(d) can be extended to an integer function of d for real (and 
complex) d. Can one say that for generic real d the volume of the unitary sphere is given by 
V{d)l 

If an unitary sphere in a space with an non-integer dimension would exist, it is likely that 
its volume would be given by eq. ([!]). Unfortunately a space with an non-integer dimension 
is not defined and therefore the question of the volume of the unitary sphere is ill-posed. 

Now we could try to reverse the question. We would like to say that a space with non- 
integer dimensions d is characterized by the property that the volume of the unitary sphere 
is given by eq. (|I|). Obviously one integral is not sufficient to identify a space and other 
integrals must be defined (for example the volume of two or more intersecting spheres). On 
the other hands this construction is likely to be empty; for example for d = —3, V(d) = —2ir 
is negative and this does not make sense if the integration measure is positive Q. 

The program may be the following; one can define geometrical objects (e.g. spheres, 
lines, planes) only giving some rules on how to compute the measure of their intersections, 
without making any reference to the coordinates of the space where they are embedded. If 
the construction is sufficiently rich, we could recover most of the geometrical results without 
having to introduce the space. Nothing forbids that some geometrical objects have a negative 
measure as far as they do no correspond anymore to real objects in a real space. 

In the following I will firstly present a tentative construction of the space with non integer 
dimensions along these lines. 



2.1 The construction of a non- integer dimensions space 

The first objects I would like to define are the functions of one (or more) vectors that belong 
to a space with non-integer dimensions. This is standard construction in physics and it has 
been discussed with full mathematical rigour by Etingof 

In the usual integer dimensional case there is no difficulty in defining a mapping from 
R n to the real numbers. If n is non-integer there is no reasonable way to define a generic 

x As we shall see later, negative dimensional spaces can be interpreted as positive dimensional spaces where 
some of the coordinates are anticommuting and the integration measure d d x is not positive definite || . 
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point x that belongs to this space, so we cannot define a function F(x) and the question of 
computing 

J d n x F(x) (2) 

does not make sense if n is non-integer. 

An alternative approach consists in giving some minimal definition of the function F(x) 
and defining / d n x F(x) as a linear functional on the appropriate space without entering in 
the definition of x. 

Let us consider a simple example of this strategy. We consider a function F(x) of the 
form 

F(x) = f(x 2 ) . (3) 
In this case for integer dimensions we have that 



d n x Fix 



(x) = nVin) / dss^fis 2 ) , (4) 
Jo 

The previous expression defines a linear functional over the functions / that is well defined 
if f(s) goes to zero sufficiently fast when s goes to infinity. The extension to non integer n 
presents no difficulty and one obtains a meromorphic function of s if the function / is C°° at 
the origine. 

The previous case may seems rather peculiar, however we can extend it to the case where 
the function F depends on the vectors x only via scalar products with other vectors (y). 
The whole construction should be invariant under the rotation group 0(n), if we rotate the 
vectors x and y. 

The simplest case are functions of two vectors and three vectors that respectively have 
the following form: 

G(x,y) = g(x 2 ,y 2 ,x-y) , (5) 

and 

H(x, y, z) = h(x 2 , y 2 , z 2 , x ■ y, x ■ z,y ■ z) . (6) 

The integral over x of the previous functions is obviously well defined for integer n and we 
will show that the final formulae are defined also for non-integer n. 
Let us consider the following expression 



d n y g{y 2 ,x-y) , (7) 
evaluated for x 2 = r' 2 . In the usual case, going to radial coordinates, one finds that 

roc rir 

d n y g(y 2 , x ■ y) = nV(n - 1) / dss^ 1 / d9 sm(6) n f(s 2 , sr cos(0)) . (8) 



The previous expression makes sense also for non-integer n. In many cases it is a meromorphic 
function of n and this allow us to define the r.h.s. of equation (|S|) also in the region where 
the integral is not convergent |4]]. 

This construction can be generalized to more complex cases; however an alternative and 
simpler construction, that gives the same results, is the following. We suppose that the 
function g has a generalized Laplace transform and that it can be written as: 

/*oo /*oo 

g(y 2 ,x-y)= da d(3 G(a,(3)exp(-ay 2 - (3x-y) . (9) 
Jo Jo 
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Using standard Gaussian integrals one finds that 



/■oo roo fiY\ n l <1 ( B 2 r 2 \ 

d n yg(y 2 ,x-y) = J o da J q d(3G{a, 0) (-) exp ( ^-J . (10) 

In the same way we may easily compute integrals of the form 

d n y g(y 2 ,y-x!,...,y-x k ) . (11) 



Generally speaking a function of k vectors depends on the k(k + 1)/2 invariants []. In other 
words / is a bona fide function of one variable, g of three variables and h of six variables. 
Only in our perverse imagination the arguments of the functions F, G and H are vectors in a 
non-integer dimensional space. The crucial point where the dimension of the space enters is 
when we attach a meaning each term of the expression integral of a function over the whole 
space. 

We notice that the integral of a positive function is not a priori positive, as can be seen 
from formula (RJ). Indeed the volume of the unit sphere is no more a positive function for 



negative dimensions. For example negative values for the integral (|Tl|) may arise if the k 
vectors Xj (i = 1, k) form an orthonormal base and the value of k is larger that the dimension 
of the space. In other words we can always decompose an n dimensional space in the product 
of an usual k dimensional space (with integer k) with a space of dimension n — k. If n — k 
is negative the corresponding measure may be no more positive definite. In this case, if we 
denote y^ — x k ■ y, we get: 

/ d n y g(y 2 , y x , y 2 , . . . y k ) = (n-k)V(n-k) J dy x dy 2 . . . dy k drr n ~ k ~ l g(r 2 + ^ y 2 , y u y 2 ---Vk) , 

i=l,k 

(12) 

where the integral is done on the region r > and it should be understood in distribution 
sense as analytic continuation from the convergence region. In the familiar case k = n we 
recover the usual formula: 

/ d n y g(y 2 , y u y 2 , . . . y k ) = / dyidy 2 . . . dy k g{ ^ y 2 , yi, y 2 . . . yu) ■ (13) 
J J i=i,k 

In a similar way we can define derivatives, Laplacians and so on. Generally speaking we 
can study vector-valued functions, under the constraint that they transform in a covariant 
way under the rotational group: e.g. 

xj(x 2 ). (14) 

In other words we study our space using the filter of rotational invariance: we consider 
only rotational invariant (or covariant) functions and we define a set of transformations (e.g. 
integrals and derivatives) that depend analytically on a parameter n; for integer values of n 
these transformations coincide with the usual ones for functions defined on a n-dimensional 
space. 



2 For integer n and sufficiently large k the k(k + l)/2 invariants are not independent quantities, however 
this fact does not concern us directly, because we are considering a generic dimensional space. 
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In general the geometrical meaning of this construction is unclear. The only case where 
there is a simple geometrical interpretation is the case of negative integer dimensions. For 
example in the case of n = —2 we have that 

J d~ 2 xF(x) = -2tt J d6d6F(6) , (15) 

where 9 is an anticommuting variable. Indeed in the case where F(x) = f(x 2 ) both integrals 
give 

df 



dx 2 



x 2 =0 



(16) 



Negative even dimensional spaces can be interpreted as spaces with anticommuting (Grass- 
mann) coordinates as it was noticed a long time ago ||. 



3 Zero-dimensional matrices 
3.1 The replica method 

The interest in zero-dimensional matrices has been triggered by the so called replica method 
[Q] for disordered systems. 

In many problems of statistical mechanics of random system we have to compute the 
quantity 

/ = lim f(N) , (17) 

where 

- f(N) = N- 1 J dfi(J) \n(Zj(N)) , (18) 

Here the partition function Zj(N) is usually a positive quantity (defined in some compli- 
cated way) and dfi(J) is a normalized measure (/ d/x(J)). For example in the Sherrington- 
Kirkpatrick model for spin glasses the partition function is given by 

Zj(N) = ^expl(3 E Ji,kWk) , (19) 
W \ i,k=l,N J 

where the sum over {a} denotes the sum over the 2 N configurations of the iV variables 
Oi that may tale only the values ±1. The measure dfi(J) is a Gaussian measure over the 
N(N + l)/2 components of the symmetric J matrix whose elements are identically distributed 
uncorrelated Gaussian variables with zero average and variance iV -1 . 
In the usual approach it is convenient to define the quantity 

- $(JV,7i) = (Nn)- 1 In (| dfi(J)(Zj(N) n ^ . (20) 

As it can be trivially seen, if the appropriate integrals are convergent, <&(N, n) is an analytic 
function of n, that is regular around n = 0, and 

f(N) = lim iMN,n) . (21) 
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This baroque construction has been introduced because in some problems for integer n one 
can perform the integration over the variables J and one can prove the following exact formula 

exp(-Nn$(N,n)) = J dQ exp(-NF n (Q)) , (22) 

where the integral is done over all n x n matrices Q a ^ {a,b = 1, n). such that its diagonal 
elements are equal to 1 (i.e. Q a , a = 1)- Using the method of the point of maximum one can 
rigorously prove that for integer n 

<j){n) = lim $(iV,n) = F n (Q*) , (23) 

N— >oo 

where the matrix Q* minimizes the function F n (Q). 

Forgetting physical motivations the mathematical problem is the following: given a func- 
tion F n (Q) defined on n x n matrices one should go through the following steps: 

1. To prove that the quantity <3>(iV, n) defined in eq. (^) can be extended to an analytic 
function of n. 

2. To prove that the limit f(N) = lim n ^o *&{N, n) exists. 

3. To prove that the quantity f(N) has a limit (/) when N goes to infinity. 

4. To give a closed expression for the quantity /. 

Of course everything depends on the form of the function F n (Q). For example in the case 
where 

F n (Q) = ^Tr(Q 2 ) (24) 
everything is trivial. One explicitly finds that 

$(iV,n) = -^ln^ (25) 

and all the issues on the dependence on n and iV of $(iV, n) can be easily settled. 

If the function $(iV, n) has also the representation eq. (|20|) , points 1 and 2 are easy, 
points 3 is more difficult, but it can be solved in some interesting cases |TB[. At the present 



time points 4 does not have a rigorous mathematical solution in the general case, 

In order to gave an idea of the problems physicists are interested, it is convenient to write 
down some of the functions that have been studied in the literature. 

F n (Q) = -Tr[ln(Q)]+/3 2 Tr(Q 2 ) 
F n (Q) = -lY[ln(Q)] + (3 2 J2(Qa, b ? (26) 

a,b 

F n (Q) = -G n (Q) + (3 2 Tr(Q 2 ) 
where in the last case (that correspond to the SK model, where we have that: 



exp(nG n (Q))= U 



a=l,n 



E 

>a=±l 



exp Qa,b(?aVb (27) 



a.b 
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The physicists have used the following approach. For integer n the matrix Q* that mini- 
mizes the function F n (Q) can be easily found and in this way one computes 4>(n) for integer 
n. Usually 4>{n) has a simple expression that can be analytically continued to n — 0. We 
would naturally guess that this analytic continuation gives the exact result. 

However in general this is not the correct procedure. The function 0(n) is defined as the 
infinite N limit of $(iV, n) and, although $(iV, n) is analytic in n, the function <p{n) may be 
a non-analytic function of n (e.g. it may have a singularity in the interval 0-1) so that the 
analytic continuation of <p(n) is not the limit for N going to infinity of $(iV, n). In some 
crucial cases is was proved that this happens (the naive approach gives inconsistencies) and 
something else should be tried |§, 

The next approach consists in trying to use the method of the point of maximum directly 
for non integer n, with all possible complications that may arise in dealing directly with 
functions defined on a non integer dimensions case M. This approach is consistent (in the 
sense that the final results are not contradictory). Its results have not yet been proved to 
be correct and goes under the name of spontaneously broken replica symmetry (the naive 
approach is usually called the replica symmetric one). 

It is a great challenge to people working in probability theory to prove (or to disprove) 



that the results of the replica broken approach are correct [12 



3.2 The replica symmetry 

We have already seem that symmetry considerations play a crucial role in defining quantities 
on non-integer dimensional spaces. It is therefore reasonable to investigate the symmetry 
properties of our problems. Moreover the presence of a symmetry group in a function F 
allow us to divide the set of its stationary points into orbits under the action of symmetry 
group. 

Let us introduce the usual permutation group S n (i.e. the replica group) and let us denote 
by ir(a) (for a = 1, n) the element where a is carried by the permutation tt. The permutation 
group (that is a subgroup of the rotational group) acts in a canonical way on the matrices 
Q, permuting both row and columns: 

Ql,b = Qn(a),n(b) ■ (28) 

In the cases we are interested the permutation group is a symmetry of the problem, i.e. 
F{Q 7T ) = F(Q). This symmetry is related to the procedure used in the construction of the 
function F it is clearly exact in the examples we have considered. 

If a function F has a symmetry, the symmetric point is a natural candidate for a minimum: 
it is automatically a stationary point, i.e. a solution of the equation 

dF . . 

■ . (29) 



dQ a 



A matrix that is left invariant by the action of the permutation group must have all off- 
diagonal elements equal: 

Q atb = q for a ^ b . (30) 
(The diagonal elements of Q are equal to one by construction.) 
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However the symmetric solution of equation fl29|) does not need to be a minimum, it may 
be a maximum, or a saddle. In order to establish the nature of the symmetric solution we 
have to compute the eigenvalues of the Hessian 

d 2 F 

%a,b;c,d = T^j J\ • (31) 

In one of the most interesting cases (the Sherrington-Kirkpatrick model for spin glasses M, 
i.e. the last of our examples) one finds that the eigenvalues of the Hessian are all positive for 
positive integer n, their analytic continuation is still positive for n > n c , but some of them 
change sign at n c where < n c < 1 0. By changing the parameter n a minimum becomes 
a saddle point and an other stationary point becomes the true minimum: this a very usual 
phenomenon (e.g. a bifurcation or spontaneous symmetry breaking) in an unusual setting. 



3.3 Breaking the replica symmetry 

If the symmetric point is not a minimum, we can look for a minimum that is not symmetric. 
In this case physicists say that the symmetry group is spontaneously broken. The subgroup 
that leaves invariant the true minimum (the little group) is the surviving (or unbroken) 
symmetry group. 

Our task is to find a new candidate for a minimum. Different attempts have been done. 
The most promising candidate ||, || is given by the matrix Q, whose off-diagonal elements 
can be written as 

Q a ,b = qi if I(a/m) = I(b/m) , 

Q a , b = g if I(a/m) ? I(b/m) , (32) 

where I(x) is the integer part of x and m is an integer that divides n. Obviously if q\ = q = q 
we recover the previous case. The symmetry group that leaves invariant this matrix is (for 
Qi Qo) the semidirect product of S n / m with the direct product of n/m copies of S m and it 
is obviously a subgroup of S n . 

The value of F(Q) evaluated at this saddle point gives a function F(q ,qi,m), that can 
be written in a relatively simple form, depending on the model (in the worst case it involves 
a few integrals). 

Now a natural question is which values of m are allowed. Originally m should be an 
integer that divides n, however, when n is no more an integer, there is apparently no reason 
for imposing that m remains an integer. It is clear that we are in a nobody land and any 
prescription is arbitrary. We will present the standard prescription, that may justified a 
posteriori by its success in giving sensible results, that are quite likely the correct ones. We 
will follow the rules: 

• Non-integer values of m are allowed as soon as n is no more a positive integer. 

• The interesting region for m is 1 < m < n for n > 1 and n < m < 1 for n < 1. 

• We look for a stationary point that satisfies the condition that all the eigenvalues of Ti 
are positive in order to find the minimum of F(Q). In the previous defined interesting 
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region (n < m < 1 for n < 1) this condition implies that F(q ,qi,m) is a maximum 
as function of go an d qi- This confusing situation is due to the fact the if the number 
of dimensions is negative, the two conditions (7i has non negative eigenvalues and 
(x, Hx) > 0) are not equivalent. 

• By uniformity we decide that the function F(q ,qi,m) should also maximized with 
respect to m. The final prescription consists in finding the maximum with respect all 
the three arguments. 

This prescription is called one step replica symmetry breaking. In certain models it 
is believed to give the exact results, i.e. in the cases for which one obtains a non-negative 
Hessian (an alternative approach to these models is sketched in the appendix). In other cases 
(e.g. in the SK model) one does not obtain a non negative value of the Hessian although 
the absolute value of the negative eigenvalue decrease be a large factor (around 10). The 
prescription we have introduced is apparently able to compute in an exact way the quantity 
/ in some models. In the other models, where it does not give a consistent result, it is a step 
toward the correct solution, that will be found in the next section. 



3.4 Full replica symmetry breaking 

As we have mentioned before the one step replica symmetry scheme does not produce con- 
sistent results in the case of the SK model. If we pay attention to the scheme that we have 
used to break the permutation group, the statements of the previous section corresponds to 
saying that in the limit n — > the unbroken subgroup of S n is the semidirect product of 
the group So with the direct product of zero copies of the group S m . In other words the 
group Sq contains itself as a subgroup and it is therefore an infinite group (it also contains 
infinitesimal elements [|Kj)- This fact may be not surprising if we ponder the fact the group 
So is obtained as analytic continuation of all the S n groups, with arbitrary n. 

One can generalize the formula eq. (|32D to the case where there are k + 1 parameter q^ 
and k parameter m^. The previous formulae correspond to the case k — 1. We present the 
explicit construction in the case k = 2. We assume the n is a multiple of m 2 and m 2 is a 
multiple of rri\. The off-diagonal elements of the matrix Q, can be written as 

Q a ,b = <?2 if I{a/m 2 ) = I(b/m 2 ) , 
Qa,b = <?i if I(a/m 2 ) ^ I(b/m 2 ) and I(a/nii) = I(b/nii) , (33) 

Q a ,b = qo if I{a/m x ) ^ /(6/mi) , 
It is remarkable that after a permutation one could write write such a matrix under the form 



Qa,b = Q( a ~ b), where the function Q(a) is a simple function of a [|TTJ. Indeed this happens 
if we set 

Q(a) = q 2 if n/m 2 divides a , 
Q(a) = q\ if n/mi divides a and n/m 2 does not divide a, (34) 

Q(a) = q if n/mi does not divide a, 

It is not clear if this translational invariant formulations, that is at the heart of the p-adic 
approach ||TTJ , has a deeper meaning (up to this moment it has been used only at the technical 
level). 
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Also in this case we first compute the function F for integer n, m x and m 2 and we 
analytically continue the result up to n=0 []. At the end of the computation we do not 
require any more that the m's are integer. The final expression is maximized as function of 
the q's and of the m's, that turn out to be in the region < ni\ < 777.2 < 1. In the case where 
the one step solution is not good, (i.e. negative Hessian) also the two steps solution is not 
good (the Hessian is still negative). 

It is clear that this construction can go on for any number of steps we want in a recursive 
way. A detailed computation shows that we always find the the Hessian is negative (for 
example in the SK model near the critical temperature the minimum eigenvalue of the Hessian 
is —A(2k + 1)~ 2 with positive A). It seems reasonable that the interesting case correspond 
to the limit k — > oo. 

In order to control better the k — > oo it is important that for finite k and in the limit 
n — > the interesting solutions of the stationary equation dFj dQ a ,b are always in the region 

= m < mi ... < iTik < 1 = nik+i . (35) 

One can thus associate to the matrix Q a piecewise constant function q(x) such that 

q(x) = qi for rrii < x < m i+ i . (36) 

If we send k to infinity, one finds by an esplicite computation that the function q(x) becomes 
a continuous function. In this way we have a canonical way to associate to a function q(x) a 
matrix in the space of zero by zero matrices. 

The introduction of the function q(x) is not so arbitrary, Indeed simple expressions can 
be obtained for the ^-invariants in the limit of n going to zero also for finite k. I present 
just two simple examples: 

lim n~ Tr(Q 2 ) = 1 — / dxq(x) 2 , 
n->o " Jo 

f 1 3 f 1 f 1 1 f 1 

lim n _1 Tr(Q 3 ) = 1 — 3 / dxq(x) 2 + - dx dyq(x)q(y) 2 + - / dxxq(x) 3 . (37) 

n->0 " Jo 2 JO Jx 2 Jo 

The last formula is related to the formula for the volume of the intersection of three p-adic 
spheres [jLlJ| . 

A more complex example is given by the quantity G(0), defined in eq. (0): a non trivial 
computation implies that 

G(0) = /(0,0), (38) 

where the function f(x, h) is defined in the strip < x < 1, — oo < h < oo and it satisfies 
the following non-linear differential equation : 

?l--l*(W +x (?l) 2 ) (39) 
dx " 2dx [dh* + \dh) ) ■ 1 ' 

the with the boundary conditions 

f(l,h) = ln(cosh(/3/i)). (40) 

3 The group properties of this solution are discussed in details in |]l4j|. 
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In all the previous cases the end of the computation the free energy becomes a functional 
of the function q(x) (F[q\) and this functional has an explicit expression that in some cases 
may be rather complex. Maximizing the functional F[q] (and depending on the temperature 
range this can be done or analytically or numerically) one finds a consistent solution of many 
models (including the SK model) such that all the infinite eigenvalues of the Hessian Q are 
non negative. Moreover this solutions satisfies all inequalities that one derive from general 
principles (e.g. positivity of the entropy). 

It obvious that all this does not make sense from a strict mathematical point of view be- 
cause the prescriptions we have used are somewhat arbitrary, although they look reasonable. 
On the other ends the replica method has been very successfully as can be seen in the next 
section. 



4 The success of the replica method 

In the case of the Sherrington Kirkpatrick model the first success of the replica method, 
especially in the case of full replica symmetry breaking, is the existence of a consistent 
solution (a result that is not evident a priori). 

It is also remarkable that in all the many infinite range models, where the replica method 
has applied, the results are in a very good agreement with the computer simulations and 
they coincide with the exact results, when available. 

It is very interesting that one can rederive most of the results of the replica method, 
without using replicas, but doing probabilistic arguments. It appears that the form of the 
matrix Q codes some information on the probability distribution of the equilibrium configu- 
rations of the system in a way that it would take too much time to explain || |{|. We only 
say that the organization of the pure states of the systems (called lumps by Talagrand |L2] 



in this context) mirrors the form of the matrix Q: e.g. the ultrametric form of the matrix 
q, that is a remnant of its p-adic origine, implies that the distances among states satisfy 
the ultrametric inequality. This alternative probabilist approach is not evident and it was 
found a few years after the establishment of the replica approach. Moreover in some models 
non-trivial rigorously results have been obtained that point toward the correctness of the 



results obtained via the replica method [12 



Most recently in a very clever paper Guerra have been able to prove that in the SK model 
(and in similar models) the following inequality si satisfied: 

f>F[q] Vq(x). (41) 

This result is important because it shows that the functional F[q], that is obtained via the 
replica method, is related to the true value (/) of the free energy. 

The same results of the replica method may be obtained in the framework of off-equilibrium 
dynamics that has been shown to be equivalent to the replica method [|i~5| , |j~6f . Recent ex- 
perimental results |T7J beautifully confirm the theoretical predictions of this approach. 

4 The Hessian becomes now a rather complex integral operator and it has both an infinite point spectrum 
and a continuous spectrum. 
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The existence of alternative routes clearly shows that the results make sense. In any case 
the replica method has a very strong heuristic value and it would be surprising if we could 
not assign a precise mathematical meaning to such an useful method. 



Appendix 

In this appendix we report an explicit computation in the case of one step replica symmetry 
breaking that may have some chance to be the starting point of a rigorous mathematical 
approach [|T8 1 . As we shall see, this result suggests that there is a correspondence from a 



class of x matrices with an infinite set ofnxn matrices for any n. 

At the technical level the idea is relatively simple. We write for negative s 

Y(-s)Z s = / dtr^expi-tZ) = / dtt~ l ~ s Y ^-Z n , (42) 

where everything has an implicit dependence on N and the bar denotes the average over J. 
If we continue the previous formula up to s —>■ we obtain and expression for ln(Z s ) 

If we can control the expression for Z n with sufficient accuracy when N — > oo, we have 
done our job. This task is not easy, because of possible cancellations due to the minus signs: 
the leading term of the sum may be not the sum of the leading terms for Z n and therefore 
we have to control Z n with a very high accuracy. 

We have already seen that the quantity Z n has an exact expression that is given by eq. 
(p2[). An interesting conjecture is that we obtain the asymptotically correct result for Z s , if 
inside eq. ( f4"2"D we make the approximation 

^«£exp(-jVF n (Q*)) , (43) 

where the sum is done over all the n x n matrices Q* that are local minima of F n (Q). This 
approximation obviously gives the asymptotically dominant result for Z n , where only the 
global minimum counts, however it is not evident that it gives the correct result when it is 
used for computing Z s . 

In the general case we do not know if conjecture is true, how to compute all the mimima 
and if this conjecture does reproduce the result of the replica approach. However in some 
models where replica symmetry is broken at one step, the appropriate computations can be 



done [0 and one finds the same result obtained with the replica method. 

There is indeed a class of models where for positive integer n all the minima of F[Q] can 
be obtained in the following way: one divides the indices in M sets of size (k = 1, M) 
where J2k m k = n. One sets 

Q a ,b = qk , (44) 

if both a and b belongs to the k th set; otherwise one has Q a ^ = 0. In some cases |1| one 
finds that if the leading contribution (that coming from the global minimum) for each n is 
used, the result coincides with the one coming form the replica symmetric approximation 
(that for a certain range of the parameters is not correct). One the other end, if one consider 
the contribution coming from all possible minima (that are subdominant for give n), at the 
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end one recovers the replica broken result where one has to find the maximum of the free 
energy for m in the region < m < 1. 

We have no idea of how the afore-mentioned conjecture can be proved and it would be 
interesting to test it in a wider context. We have mentioned it here because it hints that 
the matrices for non integer n and m used in the replica method are related to families of 
matrices that are defined for all integer n. 
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